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Introduction 


This  report  gives  an  overview  of  the  strueture  and  performanee  of  LAGER  Opties  vl.O  (Loeal 
Automated  Glider  Editing  Routine;  Optics,  version  1),  using  field  exercises  data.  EAGER  Optics 
vl.O  is  the  first  version  of  a  suite  of  algorithms  for  automated  data  processing  from  optical 
sensors  within  the  ocean  glider  quality  control  system  developed  for  the  Naval  Oceanographic 
Office  (NAVO)  Glider  Operations  Center  (GOC).  The  EAGER  Optics  vl.O  is  an  upgrade  to  the 
successful  EAGER  vl.O,  which  has  been  transitioned  to  NAVO  and  currently  operational  in 
NAVO  GOC.  The  EAGER  vl.O  methods  and  operation  are  fully  documented  in  Carnes  (2008) 
and  Carnes  and  Hogan  (2009).  EAGER  Optics  vl.O  is  a  core  part  of  EAGER  v2.0,  which  also 
includes  algorithms  updates  and  conversion  from  previous  Matlab  functions  to  EORTRAN 
subroutines.  The  complete  description  of  EAGER  2.0  will  be  found  in  a  separate  document. 
EAGER  vl.O  was  designed  and  developed  primarily  to  process  data  from  temperature  and 
conductivity  sensors  (from  which  salinity  and  density  are  derived).  It  was  developed  to  import 
raw  data  returned  from  ocean  gliders  operated  by  NAVO,  to  detect  and  flag  erroneous 
observations,  and  to  reformat  data  and  flags  into  a  format  accepted  by  the  NAVO  Real-Time 
Data  Handling  System  (RTDHS)  in  Binary  Universal  Eormat  for  the  Representation  of 
Meteorological  Data  (BUER,  details  below),  adapted  for  oceanographic  data  types.  The 
apparently  erroneous  glider  observations,  which  cannot  be  corrected  automatically,  require 
manual  intervention.  An  interactive  GUI-based  manual  utility  GUI  (MUG)  for  man-in-the-loop 
quality  control  is  used  to  change  the  quality  control  (QC)  flags  and/or  to  edit  the  observations 
using  the  controls  available  in  the  GUI,  operated  by  NAVO  GOC  personnel.  The  initial  version 
is  written  in  Matlab  (MathWorks  Inc.).  A  EORTRAN  version  of  MUG,  is  under  development, 
and  will  be  delivered  under  EAGER  v2.5.  This  validation  test  report  is  primarily  aimed  to  test 
the  automated  optics  QC  module,  which  was  also  first  developed  in  Matlab  and  subsequently 
rewritten  in  EORTRAN  77/90  to  enable  faster  processing  in  an  operational  environment. 

To  answer  the  needs  of  Mine  Warfare  (MIW),  Anti-submarine  Warfare  (ASW),  safety  of 
navigation,  monitoring  global  climate  change,  and  battlespace  environment  sensing,  optical 
sensors  have  been  fitted  and  tested  on  different  types  of  gliders,  and  have  proven  successful  in 
assessing  environmental  conditions.  They  provide  inputs  as  well  as  validation  to  system 
performance  prediction,  diver  visibilities,  optimizing  3-dimensional  ocean  optical  forecasting, 
and  benefit  higher  level  tactical  decision  making.  The  data  flow  can  be  viewed  as  a  5-step 
process,  with  the  inclusion  of  optics,  as  shown  in  Eigure  1. 

When  fully  implemented,  EAGER  v2.0  will  incorporate  the  Optics  quality  control  (QC) 
algorithm  suite  described  in  this  report.  The  optics  algorithms  will  be  used  to  process  data  from 
several  optical  sensor  types,  including  instruments  that  measure  the  beam  attenuation  coefficient 
(c),  the  total  scattering  coefficient  {b),  the  backscatter  coefficient  {bb),  the  fluorescence  returns 
and  thus  derived  concentrations  of  chlorophyll,  phycoerythrin,  and  colored  dissolved  organic 
matter  (CDOM),  downwelling  irradiance  {Ed),  phytosynthetically  available  radiation  (PAR),  and 
other  derived  optical  parameters  such  as  diver  visibility  (Table  1).  All  current  optical  sensors 
including  ECO  series  from  Wetlabs  (bb2f,  bb2slo,  fi3,  SAM,  BAM,  AUVb,  ECO-PAR)  and 
OCR504/7  from  Satlantic  are  included  in  the  current  tested  version,  with  flexibility  built  in  to 
allow  for  expansion  to  meet  future  needs.  Additionally,  Table  1  outlines  the  optical  parameter 
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typical  range,  resolution,  bandwidth,  and  channel  limits,  whieh  are  used  in  the  LAGER  Optics 
vl.O  routines  and  will  be  discussed  in  later  sections. 
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Figure  1.  Data  flow  chart  of  LAGER  Optics  vl.O.  Data  ingest;  reads  in  real-time  raw  data  at 
NAVO,  in  eonverted  netedf  format.  Auto  QC  is  the  key  step  in  LAGER  Opties  vl.O  and  the  main 
subjeet  of  this  test  report.  It  uses  flags  to  classify  different  types  of  erroneous  data  and  the 
eombined  quality  of  the  flags  determines  whether  Manual  User  GUI  (MUG)  is  needed  for  closer 
inspection.  The  resulting  data  is  then  sent  to  the  database  in  BUFR  format  required  by  RTDHS. 

To  meet  automated  data  QC,  downstream  model  requirements  and  arehive  needs,  BUFR  format 
is  adapted  from  meteorological  research  to  be  used  at  NAVO  GOC.  All  glider  sensor  outputs 
ineluding  temperature,  salinity  are  examined  and  converted  in  BUFR  data  tables.  Details  of  such 
tables  can  be  found  in  Mangin  2009  (Ref  3).  The  NRF/NAVO  optics  table  (Table  1)  has  been 
derived  and  incorporated  into  NAVO  BUFR  tables,  with  an  assigned  Class  60.  For  convenience 
of  discussion,  the  optical  variables  corresponding  to  Table  1  are  shown  in  Table  2,  adopted  from 
Ref  3  for  eonvenience. 

This  report  describes  and  validates  the  automated  LAGER  V2.0  (Optics  Vl.O)  algorithms  using 
backscatter,  fluorescence  and  irradiance  observations  from  several  Slocums  and  Seagliders 
operated  by  NAVO  and  by  NRL  during  recent  exercises  and  research  experiments.  During  the 
recent  US  Navy  Arabian  Gauntlet  2009  exereise,  EAGER  was  used  operationally  by  the  NAVO 
GOC  to  test  its  performance.  The  validation  reported  here  is  primarily  intended  to  demonstrate 
that  the  EAGER  seheme  for  deciding  whether  to  send  an  optics  profile  for  manual  QC,  based  on 
the  automated  QC  results,  is  adequate.  The  scheme  works  well  if  it  sends  profiles  to  manual  QC 
only  when  there  are  unflagged  errors  remaining  in  the  file  after  the  automated  QC  operations  are 
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applied.  A  complementary  validation  of  the  hydrographic  observation  QC  results  has  already 
been  performed  as  part  of  the  LAGER  1.0  reporting  process  and  is  described  in  Carnes  and 
Hogan^  (2009). 
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Figure  2.  LAGER  Optics  Vl.O  data  quality  check  flow  chart.  Notice  that  manual  QC  in 
FORTRAN  is  currently  in  beta  stage  and  will  be  delivered  in  LAGER  V2.5. 
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2  Testing  areas 

Data  from  8  glider  deployments  at  different  loeations  eonsisting  of  7,067  profiles  (Table  4),  from 
Biospaee  (2008),  RIMPAC  2008  exereise,  Arabian  Gauntlet  2009  exereise  are  used  for  statistieal 
analysis,  necessary  for  overall  program  parameter  settings.  Exercise  Arabian  Gauntlet  2009  and 
Exercise  RIMPAC  2008  deployments  are  used  to  validate  the  EAGER  Optics  Vl.O  AutoQC 
module.  Detailed  area  maps  can  be  found  in  subsection  7  below. 

3  Calibration 

Calibration  of  the  optics  variables  is  generally  performed  prior  to  the  optics  QC  steps.  However, 
in  cases  where  the  process  optics.f  subroutine  determines  that  calibrated  data  are  already 
available  (i.e.,  the  sensor  data  was  calibrated  during  deployment,  or  immediately  after  glider 
recovery),  that  step  is  skipped.  The  calibrate  optics. f  subroutine  called  by  process  optics.f  loads 
calibration  data  from  files  stored  in  the  EAGER/SETUP  EIEES/calibration  folder,  based  on  the 
type  of  sensor  variable  (e.g.,  backscatter,  bb,  downwelling  irradiance  ed  and  such.),  its  frequency 
band  (e.g.,  470,  532,  670,  700nm,  etc.)  and  the  date  of  the  dive  (the  calibration  files  contain 
calibration  parameter  data  pertaining  to  a  range  of  dates,  spanning  the  intervals  between 
successive  calibrations).  Eor  most  sensors,  only  a  simple  linear  calibration  formula  operating  on 
the  sensor  raw  counts  and  parameterized  by  a  scale  factor  and  offset  (typically  dark  counts)  is 
needed.  This  contrasts  with  the  calibration  of  the  hydrographic  data,  which  requires  a  larger 
number  of  parameters,  specifying  higher-order  polynomials,  for  determining  temperature  and 
conductivity,  with  the  additional  complications  that  the  salinity  is  a  function  of  both  temperature 
and  conductivity,  and  there  are  significant  response  lags  that  must  be  taken  into  account  to  obtain 
valid  profiles^  (Carnes  2008).  Euture  versions  will  apply  more  sophisticated  approaches  in 
calibration  and  correction  algorithms,  including  but  not  limiting  to  fixes  to  compensate  for 
surface  fluorescence  signals. 

4  Optics  Quality  Control  tests  and  flag  setting  at  each  depth 

Quality  control  tests  for  optics  variables  implemented  in  EAGER,  follow  certain  algorithms  and 
conventions  for  detecting  and  flagging  bad  data  that  were  developed  for  processing  the 
hydrographic  eg.,  temperature,  T,  and  salinity,  S,  data.  However,  because  of  the  multiplicity  of 
optics  variable  types  and  frequency  bands,  a  2D  variable  array  structure  is  used  in  contrast  to  the 
ID  vector  approach  used  individually  for  T  and  S.  The  reader  is  referred  to  the  EAGER  1.0 
manual'  (Carnes,  2008)  for  details  of  the  hydrographic  data  QC  implementation.  Two  of  the 
hydrographic  quality  control  tests  (algorithms)  implemented  in  LAGER  are  used  with  minor 
adaptations  for  the  Optics  QC  (OQC).  These  are  the  Global  bounds  check  and  the  Spike  test.  The 
latter  defines  a  spike  to  be  a  single  datum  departing  significantly  from  its  neighboring  values 
(single-point  spike).  It  has  the  additional  feature  of  being  able  to  detect  such  a  spike  in  the 
presence  of  a  gradient  (as  a  function  of  depth)  of  the  relevant  hydrographic  or  optical  variable. 
Since  segments  of  anomalous  optical  data  may  span  several  data  points,  a  method  to  detect  these 
cases,  in  the  form  of  a  running  standard  deviation  filter,  is  also  included  in  the  OQC.  This  filter 
finds  and  flags  data  that  depart  significantly  in  value  from  a  local  mean  computed  inside  a 
‘running’  window  (depth  interval)  that  is  moved  through  the  entire  sampled  depth  range  of  the 
profde. 
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In  addition  to  the  variable  value  range  cheek,  the  depth  of  each  sample  point  is  also  checked  and 
‘chopped’  short,  if  necessary,  to  eliminate  values  that  appear  to  lie  above,  or  too  close  to  the 
surface,  where  processes  such  as  wind  waves  and  bubbles  could  make  optical  measurements 
such  as  backscattering  and  downwelling  irradiance  either  invalid  or  difficult  to  interpret  reliably. 
If  this  depth  is  shallower  than  a  given  depth  (default  1  m),  or  is  negative,  the  variable  value  is 
flagged  accordingly.  A  constant  profile  check  is  also  done  to  determine  if  the  values  for  a 
particular  optics  variable  are  constant  throughout  the  sampled  depth  range,  which  was  proven 
critical  in  assessing  real-time  optical  sensor  performance  such  as  those  during  the  RIMPAC  2008 
exercise.  This  might  occur  for  a  variety  of  reasons,  such  as  instrument  sampling  faults,  or  sensor 
contamination,  which  suggest  the  data  are  invalid.  Details  concerning  the  implementation  of 
these  filters  are  given  in  the  following  subsections. 

In  the  present  LAGER  Optics  1 .0  version,  all  the  OQC  tests  are  independent  of  geographic 
location  and  time  of  year,  and  the  tests  are  also  independent  of  depth,  except  in  cases  where 
different  critical  test  values  are  used  in  two  different  depth  ranges.  The  universal  character  of  the 
tests  weakens  their  capability  to  detect  erroneous  anomalies  in  specific  regions.  However,  most 
of  the  parameters  used  to  configure  the  OQC  algorithms  are  stored  in  information  files,  and  may 
be  edited  to  adapt  them  manually,  if  necessary,  to  specific  regions,  or  seasons.  A  future  version 
of  LAGER  will  include  capabilities  to  adapt  these  parameters  to  specific  regions  automatically, 
based  on  a  regional  climatology  or  historical  data,  or  alternatively,  specific  regimes,  such  as 
shallow  or  deep  ocean  settings,  or  proximity  to  the  coast,  and  relationship  to  other  physical 
parameters.  This  is  particularly  important  for  optical  data,  for  which  the  dynamic  range  varies 
widely,  depending  upon  the  prevailing  oceanographic  setting  or  conditions.  In  highly  productive, 
nutrient-rich  regions  and/or  close  to  estuarine  freshwater  sources,  values  of  Chlorophyll-a  and 
Colored  Dissolved  Organic  Matter  (CDOM)  at  the  surface,  or  in  the  pycnocline,  may  reach 
highly-elevated  values,  and  may  exhibit  very  steep  property  gradients,  that  are  not  observed  in 
deeper  water. 

Eor  the  hydrographic  variables,  several  glider-specific  tests  have  been  added  to  the  OQC  to 
detect  and  flag  specific  known  types  of  bad  behavior  exhibited  by  either  specific  brands  of 
gliders  or  by  all  types  of  gliders.  In  most  cases,  the  glider-specific  tests  are  functions  of  the 
vertical  velocity  of  the  glider  which  is  employed  as  a  substitute  for  the  more-difficult-to- 
determine  total  speed  of  the  glider  through  the  water.  There  are  some  early  indications  that 
similar  provisions  will  need  to  be  made  for  the  optics  data.  However,  insufficient  data  have  been 
accumulated  to  clearly  establish  the  need  for  such  corrections,  and  to  reliably  specify  a  remedial 
strategy. 

Optical  error  flags  are  initially  set  to  zero  at  each  depth.  Tests  are  performed  on  observed 
temperature  and  salinity  values  prior  to  any  tests  performed  on  optical  values.  If  a  test  is  failed  at 
a  specified  depth,  then  the  corresponding  flag  is  set  to  the  failure  flag  value  for  that  test  at  that 
depth.  The  sequence  of  tests  is  described  next. 

4. 1  Global  bounds  check 

The  algorithm  for  this  check  is  implemented  in  the  EORTRAN  subroutine,  process  optics.f 
V_flag(i)  =  2  if  V(i)  <  Vmin  or  V(i)  >  Vmax,  where  V(i)  denotes  the  value  of  a  particular  optics 
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variable  (e.g.,  backscatter  at  660  nm,  or  bb660),  at  depth  Z(i),  i  is  the  depth  index,  and  Vmin  and 
Vmax  denote  the  extreme  values  representing  the  aeeeptable  natural  range  of  values  of  that 
variable.  The  values  of  parameters  Vmin  and  Vmax  are  defined  in  the  /LAGER/SETUP_FILES 
optics_variables_info.dat  file  in  the  LAGER/SETEIP  EILES  folder  (shown  in  Appendix  A), 
adapted  from  Table  1  and  2.  Note  that  V_fiag=l  is  not  used. 

4.2  Surface  Depth  test 

The  algorithm  for  this  check  is  implemented  in  the  FORTRAN  subroutine,  chop  optics.f,  called 
by  process  optics.f  The  flag  value  is  set  to  7  (V  fiag(i)  =  7),  if  depth(i)  <  Dmax,  where  Dmax  is 
the  maximum  chop  depth,  maxchopdepth,  defined  in  the  process  optics.f  subroutine.  This  will 
be  moved  to  the  optics_variables_info.dat  file,  in  a  subsequent  version. 


4.3  Spike  test 

The  algorithm  for  this  check  is  implemented  in  the  FORTRAN  subroutine,  de  spike  time  called 
by  process_optics.f  The  flag  value  is  set  to  4(V_fiag(i)  =  4),  if  |V(i)-(V(i+l)+V(i-l))/2|  - 
|V(i+l)-V(i-l)|  >  K,  where  K  =  Ks  (Z(i)  <  500  m)  or  K  =  Kd,  (Z(i)  >=  500  m).  Ks  and  Kd  are  the 
minimum  deviations  that  could  represent  a  spike  in  the  depth  ranges  0-500m  or  >=  500m, 
respectively.  Notice  that  the  spike  test  utilizes  the  optical  parameter  table  (Table  1),  particularly 
the  resolution  of  the  instruments  for  this  task.  Changes  can  be  made  through  values  defined  in 
EAGER/SETUP  FILES/optics  variables  info.dat  (Appendix  A). 

4.4  Running  Standard  Deviation  test 

The  algorithm  for  this  check  is  implemented  in  the  FORTRAN  subroutine,  stdespikesc.f  called 
by  process_optics.f  V_flag(i)  =  5  if  |(V(i)-Vmu(i-w/2:i+w/2)|  >  F  where  F  =  stdevfactor*Vsd(i:- 
w/2;i+w/2),  w  is  the  width  of  the  moving  window,  Vmu  and  Vsd  are  the  local  running  mean  and 
standard  deviation,  respectively,  of  the  data  falling  within  this  window,  and  stdevfactor  is  a 
specified  fraction  or  multiple  of  Vsd  .  If  however,  the  absolute  difference  |(V(i)-Vmu(i- 
w/2:i+w/2)|  <  abslimfactor,  where  abslimfactor  is  a  specified  multiple  of  the  sensor  resolution, 
then  the  datum  is  not  flagged.  The  rationale  for  this  is  that  deviations  that  do  not  depart 
significantly  from  the  running  mean  may  be  retained  as  they  could  be  accounted  for  simply  by 
normal  instrument  noise  fluctuations.  They  will,  in  any  case,  be  eliminated  or  reduced  in 
subsequent  smoothing  operations.  The  values  of  stdevfactor  and  abslimfactor  are  specified  in  the 
optics_variables_info.dat  file  which  defines  QC  parameters  for  each  variable  type. 

This  filter  is  adaptive  in  that  it  detects  values  that  are  anomalous  to  the  extent  that  they  deviate 
significantly  from  the  local  mean  value,  with  ‘significantly’  referring  to  a  large  departure,  or 
outlier,  relative  to  the  local  variability  of  the  variable,  or  normal  instrument  noise.  Because  the 
center  of  the  running  window  moves  through  successive  depth  indices,  i,  each  datum  is  assessed 
in  the  context  of  its  spanning  window.  It  is  thus  possible  for  adjoining  data  values  to  be 
considered  anomalous  (outliers)  and  to  be  flagged  accordingly. 


6 


A  special  provision  is  made  in  this  test  to  ensure  anomalous  values  are  detected  near  the  very  top 
and  bottom  of  the  profile.  This  is  done  by  using  a  fixed  window,  within  which  all  data  are  tested, 
at  the  top  and  bottom  of  the  profile,  while  a  running  window,  centered  on  each  datum  in 
succession  is  used  at  intermediate  depths.  As  a  consequence  the  test  is  biased  slightly  deep  and 
shallow,  respectively,  near  the  top  and  bottom  of  the  profile,  but  this  is  considered  a  minor 
tradeoff  considering  the  importance  of  testing  all  data  values  near  these  depth  extremes. 

4.5  Constant  Profile  test 

The  algorithm  for  this  check  is  implemented  in  the  FORTRAN  subroutine,  const_profile.f,  called 
by  process  optics.f  V  flag(i)  =  6  if  max(V(i))-min(V(i)  <  Vres,  where  Vres  is  the  resolution 
(minimum  digitized  sample  increment)  of  the  given  sensor.  The  resolution  is  defined  in  the 
optics_variables_info.dat  file  (Appendix  A).  This  is  a  critical  test  to  detect  malfunctions  of  newly 
developed  optical  sensors  and/or  deployment  errors. 


Figure  3.  Example  to  show  profiles  in  question,  using  Exercise  data  from  RIMPAC  2008.  The 
current  constant  check  alerts  the  system  and  sends  the  profile  in  question  to  manual  QC  for 
inspection.  Density  profile  is  shown  in  blue. 
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4.6  Backscatter  Buddy  check 


In  the  algorithm  process  optics.f,  there  is  a  eheek  to  see  if  all  baekseatter  bands  simultaneously 
exhibit  flag  4  (single  point  spike)  or  flag  5  (standard  deviation  anomaly)  at  the  same  depth  level. 
If  that  is  the  ease,  the  error  flags  are  turned  off  The  rationale  for  this  is  that  a  spike  or  anomaly 
oeeurring  in  all  the  bands,  is  unlikely  to  be  an  isolated  instrument  or  sensor  sampling  error,  but  is 
likely  to  be  due  to  an  aetual  observable  ehange  in  water  mass  eharaeteristies. 


4.7  Interpolation  over  spikes 

The  algorithm  interpdepth.f  ealled  by  proeess  opties.f  optionally  interpolates  linearly  over  all 
missing  data  (gaps)  in  the  profde,  as  well  as  replaeing  erroneous  (flagged)  values  with 
interpolated  values.  The  interpolation  is  linear,  with  respeet  to  depth.  The  option  to  interpolate 
over  missing  data  gaps  is  presently  turned  off,  to  limit  the  risks  involved  in  extrapolating  data 
over  large  depth  ranges,  at  points  removed  from  aetual  observations. 


4.8  Convert  beta  to  bb 

Although  not  shown  in  the  flow  ehart  (Fig.  2),  the  algorithm  beta2bb.f  is  ealled  from 
proeess  opties.f  to  eonvert  the  original  and  interpolated  baekseatter  sensor  data  from  beta, 
whieh  is  observed  at  a  partieular  baekseattering  angle  (114“  eurrently),  to  bb,  whieh  is  a  measure 
of  baekseatter  integrated  over  the  baek  seattering  angles  (90  to  180)  observed  by  the  sensor.  This 
derivation  is  presently  done  after  the  opties  QC  algorithms  are  applied.  In  a  future  version, 
eonsideration  will  be  given  to  doing  this  eonversion  prior  to  QC,  to  avoid  spurious  values 
appearing  due  to  errors  in  ealibration  offset. 


4.9  Smooth  optics  variables 

After  the  QC  algorithms  are  applied  and  interpolation  has  been  performed  to  eorreet  flagged 
values,  a  smoothing  operation  is  also  applied.  The  operation  is  performed  by  the  smoothdist.f 
subroutine,  ealled  by  proeess  opties.f  The  resulting  smoothed  variable  is  output  separately  from 
the  unsmoothed  interpolated  data.  The  algorithm  smooths  by  averaging  the  data  within  a  moving 
window  spanning  a  speeified  depth  range.  The  window  size,  whieh  is  eurrently  applied  to  all 
opties  sensor  data  is  set  to  5  m.  This  parameter  value  was  seleeted,  based  on  speetral  analysis  of 
the  baekseatter  data,  whieh  shows  that  it  eoineides  with  its  spatial  deeorrelation  seale.  The 
deeorrelation  seale  represents  the  ‘eutoff  wave  length  where  the  ‘red’  speetrum  of  the 
baekseatter  signal  levels  out  at  a  low  value  presumed  to  eoineide  with  the  instruments  noise 
floor.  In  a  future  version  of  the  eode,  this  parameter  will  be  made  aeeessible  through  the 
opties_variables_info.dat  file  (Appendix  A),  so  that  its  value  ean  be  eustomized,  if  neeessary,  for 
eaeh  individual  sensor  type. 
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5  Overall  optics  quality  flags 

As  discussed  previously,  the  following  flags  are  used  for  eurrent  opties  QC  modules: 


Flag  name 


assigned  value 


Global  bounds  eheek  2 

Bb  buddy  eheek  4 

Running  STD  eheek  5 

Constant  value  eheek  6 

Surfaee  depth  eheek  7 


The  LAGER  Opties  1.0  automated  QC  sends  a  multi-sensor  multiple-band  opties  file  for  further 
evaluation  in  the  manual  QC  step  if,  for  any  of  the  individual  profiles  in  the  file,  the  number  of 
profde  observations  flagged  as  bad  is  greater  than  a  given  percentage  (default  15%,  for  the 
eurrent  validation  test,  this  value  is  set  to  10%)  of  the  total  number  of  observations  of  that 
profile.  In  future  versions,  additional  eriteria  will  be  used  to  determine  if  profiles  are  bad  enough 
to  be  deleted  or  require  manual  editing  on  the  basis  of  the  temperature  and  salinity  or  opties 
variable  values.  (See  Ref  1  for  details  eoneeming  the  eorresponding  eriteria  for  the  hydrographic 
data).  Hence,  a  file  may  be  sent  to  manual  QC  because  either  or  both  of  temperature  or  salinity 
require  manual  editing,  or  beeause  at  least  one  of  the  optics  profiles  requires  manual  editing.  In 
addition,  an  individual  T,  S  or  opties  variable  profile  may  be  flagged  for  deletion  if  the  number 
of  bad  data  or  gaps  are  exeessive.  A  needs  manual  editing  flag  is  assigned  in  eaeh  file  to  eaeh 
ascending  or  descending  profde  (eaeh  value  summarizes  quality  of  the  eorresponding  opties 
variable  profile). 

needs  manual  editing  flag  =  0  indieates  that  manual  editing  is  not  required. 

needs  manual  editing  flag  =  1  indieates  that  the  file  eontaining  this  profile  should  be  sent  to  the 

manual  editor  for  further  examination  even  if  the  needs  manual  editing  flag  for  the  other  profile 

(if  it  exists)  in  the  file  is  set  to  0  (doesn't  need  manual  editing). 

needsmanualeditingflag  =  1  if: 

1 .  If  a  given  pereentage  (default  15%)  of  the  depths  have  V  flag  >  0,  but  not 
ineluding  V  flag  =  7. 

2.  for  Sloeum  profiles  only,  at  least  one  depth  gap  between  eonsecutive  observations 
was  >=  10  m  or  two  gaps  >=  6  m. 

6  Error  Flag  Statistics 


The  goals  of  LAGER  inelude  automatieally  deteeting  and  flagging  erroneous  observations  prior 
to  being  passed  to  the  RTDHS,  and  identifying  apparently  erroneous  glider  observations  that 
cannot  be  eorrected  automatieally,  and  thus  require  manual  intervention.  It  is  thus  desirable,  to 
automatieally  deteet  and  eorreet  as  many  of  the  errors  as  possible,  while  minimizing  deteetion 
failures  (failure  to  deteet  bad  data)  and  false  alarms  (ineorreetly  rejeeting  good  data).  Sinee 
manual  intervention  is  labor  intensive,  we  also  want  to  minimize  the  number  of  data  being  sent 
(to  MUG)  for  Manual  QC. 
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There  is,  however,  inevitably,  a  degree  of  subjeetive  judgement  involved  in  determining  the 
appropriate  parameter  to  use  in  the  AutoQC  algorithms  deseribed  above,  as  well  as  in  using 
Manual  QC  operations  to  determine,  ‘by  eye’,  whether  an  error  has  oeeurred.  This  subjeetivity 
ean  only  be  redueed  by  using  ‘expert  judgement’  to  determine  the  parameters  and  to  assess  the 
performance  of  the  QC  operations.  Hence,  that  is  the  approach  we  have  adopted  in  this  report. 

An  independent  assessment  of  physical  T  and  S  variables  is  possible  by  comparison  with 
climatological  data  bases'  (e.g.,  GDEM,  see  Carnes,  2008).  However,  such  climatologies  are  not 
yet  available  for  bio-optical  data  (backscatter,  attenuation,  fluorescence  etc).  Furthermore,  the 
bio-optical  data  exhibit  a  wider  dynamic  range,  appear  in  general  to  be  more  ‘noisy’,  and  are 
subject  to  greater  environmental  variability  than  is  the  case  for  T  and  S.  This  makes  the  optical 
variables  both  prone  to  error,  while  being  subject  to  fewer  ‘realistic’  constraints.  Determining 
whether  a  particular  profile  requires  manual  intervention  is  thus  to  some  degree  an  arbitrary 
process.  However,  as  long  as  the  ‘arbitration’  is  carried  out  by  an  expert  knowledgeable  in 
measurement  instrumentation,  and  experienced  in  acquiring  and  editing  bio-optical  data,  the 
results  can  be  expected  to  be  satisfactory. 

In  this  first  version  of  the  OQC,  we  use  a  test  based  on  a  simple  bad  data  threshold  to  determine 
if  a  profile  requires  manual  QC.  This  test  is  based  on  the  following  rationale:  If  the  percentage  of 
data  that  was  detected  and  corrected  by  the  AutoQC  process  exceeds  a  certain  threshold  (default 
10%  for  current  validation  test),  it  is  likely  there  are  other  problems  or  undetected  errors  with  the 
data  set,  that  may  require  manual  correction.  Similarly,  in  this  case,  a  manual  check  is  considered 
desirable,  to  confirm  that  the  automatic  detection  and  correction  of  a  relatively  large  number  of 
errors  was  successful. 

In  the  next  section,  we  discuss  the  validation  of  the  Automatic  and  Manual  QC  operations,  based 
on  a  range  of  data  from  various  missions  and  geographical  zones,  using  an  ‘expert  judgement’ 
approach.  In  the  remainder  of  this  section,  we  present  and  discuss  statistics  that  show  what 
percentage  of  the  data  fdes  are  sent  to  manual  QC  for  a  given  bad  data  threshold,  as  applied  to 
each  profile  and  variable  type  appearing  in  the  file.  We  also  give  statistics  indicating  the  reasons 
why  the  various  profdes  were  flagged  as  bad  (based  on  the  QC  algorithms  applied).  These 
statistics  are  applied  to  data  from  various  geographic  regions  and  they  are  accumulated  over  all 
profiles  and  all  files  obtained  from  a  specific  glider  within  the  specified  region  and  sampling 
period.  For  practical  reasons,  the  sampling  period  is  defined  to  be  a  calendar  month  or  less. 

We  note  here  that,  in  the  terminology  of  FAGFR,  a  ‘file’  contains  all  the  data  from  a  single 
‘dive’,  which  in  general  might  include  several  descending  and  ascending  profiles,  while  a  single 
profile  generally  contains  data  obtained  from  multiple  sensors  and  wavebands  (eg.  bb470,  bb532, 
fl2  cdom  &  etc.).  While  the  existing  seagliders  typically  generate  files  with  multiple  profiles,  the 
currently  available  slocums  generate  files  with  only  one,  or  a  pair  of  (ascending  and  descending) 
profdes.  In  the  existing  version  1.0  of  FAGFR,  however,  all  files,  regardless  of  glider  type,  are 
constructed  so  they  comprise  a  single  profile,  or  at  most,  a  pair  of  ascending  and  descending 
profiles.  In  a  future  version  of  FAGFR,  this  constraint  will  be  relaxed,  so  that  multiple  profiles 
(or  profile  pairs),  comprising  a  complete  ‘dive’  will  appear  in  a  single  file.  For  the  purpose  of  our 
statistical  analysis,  therefore,  a  file  generally  contains  1  or  2  profiles. 
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Results  of  the  statistieal  analysis  of  the  pereentage  number  of  files  that  are  sent  to  Manual  QC 
are  shown  in  Table  4.  This  Manual  QC  pereentage  is  based  on  the  pereentage  of  bad  data 
appearing  in  eaeh  of  the  profiles,  eonsidering  all  the  variables/wavebands  in  eaeh  profile 
appearing  in  the  file.  If  one  or  more  of  the  individual  variable/wave  band  reeords  has  more  than 
the  speeified  pereentage  threshold  of  bad  data,  then  the  eorresponding  reeords  are  marked  for 
Manual  QC,  and  the  entire  file  (dive)  is  sent  to  Manual  QC  (Mug).  Under  operational  eonditions, 
files  for  whieh  none  of  the  individual  variable/band  reeords  exeeed  the  threshold  would  bypass 
Manual  QC  and  pass  direetly  to  the  BUFR  file  ereation  proeess  (see  Carnes,  2008  for  details). 

For  the  purpose  of  the  validation  tests  diseussed  in  the  next  seetion,  however,  all  the  files  were 
sent  to  manual  QC,  and  eaeh  erroneous  profile  was  eheeked  to  determine  if  it  was,  indeed, 
flagged  for  Manual  QC  or  not. 

In  Table  5,  statisties  for  partieular  areas  and  ealendar  months  (eol  1)  and  speeifie  gliders  (eol  2) 
are  shown  along  with  the  number  of  files  eontaining  valid  opties  data  (eol  3).  A  series  of  tests 
was  performed  using  these  data  sets,  whieh  eounted  the  number  of  errors  of  any  type  used  to  flag 
the  data  appearing  at  eaeh  depth  level,  and  oeeurring  in  eaeh  variable/wave  band  in  the  file.  The 
speeified  pereentage  threshold  levels  of  bad  data  (eol  4)  was  ehosen,  not  neeessarily 
exhaustively,  from  the  range  1,  2,  5,  10,  15,  20,  25,  30  %.  The  eorresponding  pereentage  of  files 
in  the  data  set  that  would  be  sent  to  Manual  QC,  given  this  threshold  is  then  shown  (eol  5).  A 
Manual  QC  (%)  value  of  0.00  indieates  that  less  than  0.01  %  of  data  would  have  gone  to  Manual 
QC  -  a  negligible  number).  For  most  data  sets,  the  pereentage  going  to  manual  QC  is  small  when 
the  threshold  is  high  (high  toleranee  for  bad  data),  and  it  inereases  monotonieally  (but  not 
neeessarily  gradually),  as  the  threshold  is  redueed  (low  toleranee  for  bad  data).  In  exeeptional 
oases,  however,  the  Manual  QC  (%)  starts  at  a  relatively  high  level  and  inereases  as  the  threshold 
is  redueed.  The  exoeptions  apply  exolusively  to  Slooums  and,  as  will  be  discussed  below,  occur 
only  in  cases  where  a  large  number  of  profiles  are  flagged  as  being  uniform  in  value  (and  hence 
highly  likely  to  be  erroneous).  Ignoring  these  exceptions  the  Manual  QC  (%)  values  for  a 
specified  bad  data  threshold  of  15  %  are  generally  smaller  than  7%,  (with  one  exception  of 
20%)  for  Slocums  (slxxx),  and  1%  for  Seagliders  (sgxxx).  The  Seagliders  appear  to  be  generally 
less  noisy  (fewer  spikes  and  other  anomalies,  see  next  section)  and  subject  to  errors  than  the 
slocums.  However,  this  might  partly  be  a  consequence  of  sampling  bias,  since  the  Seagliders  are 
often  deployed  in  deeper  waters  (subject  to  less  variability),  further  from  the  coast  (lower  values 
and  weaker  gradients),  and  sample  over  a  wider  depth  range  (tendency  for  most  variables  to  be 
small  and  uniform  at  depth)  than  is  the  case  for  the  Slocums. 

Table  5  also  gives  an  indication  of  the  likelihood  (percentage)  of  individual  flag  values  appearing 
in  data  set,  with  the  flag  counts  (expressed  as  a  percentage  of  the  total  number  of  data) 
accumulated  over  all  variables/bands  and  files.  The  flag  types  listed  are  2  for  data  out  of  range,  4 
for  a  single  point  spike,  5  for  a  running  standard  deviation  anomaly,  6  for  a  constant  profile  and  7 
for  data  that  appears  within  1  m  of  (or  above!)  the  surface.  In  the  case  of  flag  6,  if  this  occurs  at 
all  (and  this  is  the  case  only  for  slocums),  it  is  either  small  (indicating  a  small  number  of  constant 
profiles)  or  very  large  (indicating  a  large  number  of  constant  profiles).  When  large,  this  number 
tends  to  be  exceptionally  large,  because  only  cases  where  the  entire  profile  is  uniform  are 
detected  (as  distinct  from  shorter  segments  of  a  profile,  that  otherwise  contains  good  data),  and 
because,  in  that  case,  all  data  points  in  the  profile  are  flagged  as  bad.  During  AutoQC,  when  it  is 
possible  to  interpolate  over  bad  (flagged)  data,  the  number  30  is  added  to  the  error  flag. 
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However,  for  the  purpose  of  this  analysis,  flagged  data  corresponding  to  both  interpolated  and 
uninterpolated  data  were  combined  (i.e.,  flag  values  of  32,  34  or  35,  where  counted  along  with 
flag  values  of  2,  4  and  5,  and  the  resulting  sums  were  given  the  labels  2,  4  or  5,  respectively).  In 
the  case  of  flag  types  6  and  7,  interpolation  is  never  done,  and  corresponding  data  values  are 
instead  flagged  as  missing  (and  given  special  values),  so  summation  was  not  necessary,  in  that 
case.  Note  that  data  with  flag  7  is  not  necessarily  erroneous  (it  could  have  been  correctly 
measured),  but  is  excluded  because  of  the  likely  complications  for  optical  data  interpretation  in 
the  near  surface  (surface  water  wave  and  bubble  affected  near  surface  zone). 

For  Slocums,  as  a  general  rule,  flag  values  of  2  (range  check)  appear  in  less  than  0.3  %  of  the 
data,  while  values  of  4  (single-point  spike)  can  be  as  high  as  5.3  %.  Standard  deviation 
anomalies  are  slightly  less  frequent  by  a  factor  or  2  or  more,  with  maximum  values  of  about 
2.4%.  As  indicated  above  flag  6  values  (uniform  profdes)  are  usually  small  (less  than  2.3  %),  but 
may  be  exceptionally  large  from  13  to  67  %  for  certain  regions  and  gliders.  All  of  the  slocums 
that  show  large  numbers  of  uniform  profiles  are  those  owned  either  by  NRL  (sl077,  sl079,  and 
sl083),  and  deployed  during  BIOSPACE_200806,  or  by  Rutgers  University  (sl016),  deployed 
during  Hawaii_200807  (RIMPAC  2008),  and  the  problem  might  be  accounted  for  either  by 
operational  difficulties  or  known  problems  with  internal  data  communications  that  occurred  in 
the  earlier  slocum  gliders. 

For  Seagliders,  generally,  the  percentage  of  flagged  values  is  significantly  smaller  than  is  the 
case  for  Slocums.  Flag  2  (range  checks)  occur  less  than  0.1%  of  the  time.  Flag  4  (single-point 
spikes)  occur  generally  less  than  0.5%  of  the  time  (exception  -'2.0%  for  sgl32  in  Wash_200705). 
Flag  5  (standard  deviation  anomaly  check)  occurs  typically  less  than  1%  of  the  time.  Finally,  as 
noted  above.  Flag  6  (uniform  profiles)  is  never  significant  for  Seagliders  and  Flag  7  (near  surface 
check)  occurs  less  than  2.0  %  of  the  time,  the  latter  being  similar  to  Slocums. 

It  should  be  appreciated  that  the  statistics  reported  here  are  quite  strongly  dependent  upon  the 
adoption  of  particular  values  for  the  key  parameters  employed  in  the  QC  algorithms.  As  a 
general  rule,  these  parameter  values  are  specified  in  the  optics_variables_info.dat  file  (Appendix 
A).  However,  in  some  exceptional  cases  the  values  are  directly  coded  in  the  present  version, 
while  in  a  future  version  all  of  the  significant  parameters  will  be  accessible  through  the 
information  file.  To  aid  interpretation  of  the  above  statistics,  the  principles  used  in  setting  some 
of  these  parameters  are  given  here.  A  key  parameter  value  from  which  others  are  derived  is  the 
instrument  resolution,  which  generally  accords  with  the  smallest  increment  in  value  which  can  be 
represented  in  the  digital  data  record  (This  shows  up  in  the  MUG  plots,  as  a  distinct  quantum 
jump  in  value  or  locally  uniform  values,  in  weak  gradient  areas).  It  would  be  possible,  in 
principle,  for  this  value  to  be  customized  for  particular  instruments  or  sensor  units  (via  specific 
instrument  calibration  records),  but  the  value  adopted  here,  and  in  LAGER  Optics  1.0,  is  that 
which  is  considered  typical  of  the  relevant  instrument  sensor  type  and  wavelengths.  Erom  this 
quantity  the  threshold  spike  size,  or  ‘spike’,  utilized  in  the  test  for  flag  4,  is  specified  as  a 
multiplicative  factor  of  10  larger  than  the  resolution.  The  ‘abslimfactor’  utilized  in  the  test  for 
flag  5,  is  specified  as  a  multiplicative  factor  of  3  larger  than  the  resolution.  The  ‘stdevfactor’ 
which  is  a  non-dimensional  quantity  determining  the  number  of  local  standard  deviations  that 
define  a  deviation  anomaly  for  flag  5  is  also  specified  in  the  information  file.  These  parameters 
and  others,  act  together  with  the  logic  of  the  various  QC  algorithms  to  determined  the  outcome 
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of  the  particular  QC  checks,  and  directly  affect  the  statistics  reported  here.  The  implication  is 
that  after  implementation  of  the  latest  version,  the  QC  performance  can  be  fine-tuned,  as  more 
experience  is  gained  and  larger  historical  data  sets  become  available  to  facilitate  performance 
improvements. 


7  Validation  Using  Optical  Variable  Transects  and 
Anomalous  Profiles 

The  above  sections  explain  the  LAGER  Optics  Vl.O  automated  QC  algorithms  and  statistics 
under  different  parameter  settings.  The  process  is  adjustable  using  the  parameter  table  described 
in  optics_variables_info.dat  under  /LAGER/SETUP  FILES/  directory.  It  can  be  seen  that  for  the 
7,067  profdes  tested,  statistics  show  that  about  30.5%  of  the  profdes  will  be  sent  to  GOC 
operators  for  manual  inspection,  when  a  15%  threshold  is  set  for  the  optical  data  points  in 
question  of  each  individual  profile.  The  performance  of  the  algorithm  is  further  validated  by 
visual  inspection  of  522  profiles  from  slocum  glider  sll  13  during  Arabian  Gauntlet  2009 
Exercise  (August  14-24,  2009).  The  area  is  shown  in  Figure  4,  with  the  glider  track  near  shore 
from  northwest  to  southeast,  marked  in  red. 

Continuous  profdes  from  deployments  ranging  from  August  18  to  August  19  are  used.  Figure  5 
shows  the  time  series  of  all  optical  profiles  from  backscattering  470nm  during  the  entire 
exercise.  The  selection  process  was  based  on  the  quality  of  the  optics  sensor,  the  system 
performance  reliability  including  data  collection  and  transmission,  as  well  as  statistical  analysis 
based  on  number  of  profiles  in  questions.  The  intentional  bias  is  toward  more  troublesome 
profiles,  to  test  the  robustness  of  the  algorithm.  It  is  natural  to  see  more  profiles  being  sent  to 
manual  QC  queue  for  this  reason.  Also,  to  make  sure  that  no  bad  data  slip  through  the  process,  a 
tighter  than  usual  criteria  setting  is  applied,  which  means  once  10%  of  the  data  points  in  a  profile 
is  in  question,  due  to  various  flag  settings  discussed  above,  this  profile  will  be  sent  for  manual 
inspection,  as  compared  to  default  15%  value  set.  Out  of  522  profiles  inspected  visually  by  the 
authors,  305  were  sent  to  manual  QC  for  inspection.  217  were  sent  to  RTDHS  directly.  There  are 
2  profiles  out  these  pass-thru  that  are  mistakenly  marked,  which  translate  into  0.38%  error. 

The  first  one  is  from  the  profile  number  sll  13  0818  1 151  0000401,  as  shown  in  Figure  6.  The 
questionable  behavior  is  apparent  as  the  sparseness  of  the  data  should  present  the  profile  for 
manual  inspection.  This  is  likely  due  to  the  complication  from  both  surface  check  and  lack  of 
sparseness  check.  This  can  be  improved  in  the  EAGER  Optics  V2.0. 

The  second  profile  in  question  is  shown  in  Figure  7.  This  is  another  profile  that  was  not  flagged 
to  be  sent  to  manual  QC.  Profile  is  from  Slocum  glider  sll  13  on  August  18,  2009  with  the  profile 
number  sll  13_0818_1459_0000512.  The  questionable  behavior  can  be  observed  for  data  in  the 
middle  of  the  water  column  where  single  point  spikes  were  flagged.  Such  flags  are  within  the 
criterion  stated  in  previous  sections.  Further  investigation  is  necessary  and  a  fix  will  be  in  place 
for  V2.0. 
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Chlorophyll  Concentration,  OC3  Algorithm 
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0.054 
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ch)_oc3  Code  7330/Ocean  Sciences 

Persian  Gulf  250  (MODIS-AQUA-PM)  Naval  Research  Laboratory 

Version  8  Stennis  Space  Center.  MS 


Figure  4.  Arabian  Gauntlet  2009  Exereise  area  shown  using  MODIS  image  derived  ehlorophyll 
produet  from  NRL  Oeean  Seienees  braneh  (7330).  Red  line  depiets  glider  deployment  traeks, 
from  where  validation  test  data  were  used. 
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GL_PersianGulf_200908  sl113  Backscatter  470  starting  on  2009 
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Figure  5.  Time  series  of  optieal  baekseattering  at  470nm  (bb470)  from  sloeum  sll  13  during  Arabian  Gauntlet  2009  Exereise.  The 
validation  test  performed  here  uses  data  from  August  18-19  period,  whieh  is  a  good  representation  of  overall  data  strueture,  although  it 
is  biased  towards  more  questionable  profiles.  It  has  less  data  gaps  due  to  data  transmission  problems  found  during  later  stages  of  the 
exereise. 
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Figure  6.  Sample  profile  that  was  not  fiagged  to  be  sent  to  manual  QC.  Profile  is  from  Sloeum  glider  sll  13  during  Arabian 
Gauntlet  2009  on  August  18,  2009.  The  profile  number  is  sll  13  0818  1 151  0000401.  The  questionable  behavior  is  apparent 
as  the  sparseness  of  the  data  should  present  the  data  for  manual  inspeetion.  This  is  likely  due  to  the  complieation  from  both 
surfaee  eheek  and  lack  of  sparseness  check.  Notice  individual  flag  values  are  shown  both  in  drop-down  box  on  left  panel,  as 
well  as  in  the  figure  window  on  the  right,  marked  by  green  lines. 
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Figure  7.  Sample  profile  that  was  not  fiagged  to  be  sent  to  manual  QC.  Profile  is  from  Slocum  glider  sll  13  during  Arabian 
Gauntlet  2009  on  August  18,  2009.  The  profile  number  is  sll  13_0818_1459_0000512.  The  questionable  behavior  can  be 
observed  for  data  in  the  middle  of  the  water  column  where  single  point  spikes  were  flagged.  Such  flags  are  within  the  criterion 
stated  in  previous  sections.  Further  investigation  is  necessary.  Notice  individual  flag  values  are  shown  both  in  drop-down  box 
on  left  panel,  as  well  as  in  the  figure  window  on  the  right,  marked  by  green  lines. 
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Figure  8.  Sample  profile  that  was  not  fiagged  to  be  sent  to  manual  QC.  Profile  is  from  Slocum  glider  sll  13  during  Arabian 
Gauntlet  Exercise  on  August  18,  2009.  The  profile  number  is  sll  13_0818_1358_0000504.  This  profde  has  been  flagged  to  be 
manually  examined,  although  it  does  not  show  any  apparent  errors,  after  auto  QC  process.  This  is  likely  due  to  the  tight  criteria 
currently  in  place  for  this  validation  testing. 
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Over  all,  more  profiles  were  sent  to  manual  QC  than  default  values  indicated,  due  to  the  selection 
of  certain  particular  segments  (more  troublesome  profiles).  A  majority  of  these  profiles  were  sent 
due  to  tight  spike  check  thresholds.  One  example  can  be  found  in  Figure  8.  This  profile  has  been 
flagged  to  be  manually  examined,  although  it  does  not  show  any  apparent  errors,  after  auto  QC 
process.  This  is  likely  due  to  the  tight  criteria  currently  in  place  which  should  be  relaxed  in  order 
to  lessen  the  burdens  on  the  operators. 

Conclusions 

Through  the  testing  of  7,067  NAVO  collected  glider  optics  data  profiles  from  various  Slocum 
gliders  deployed  over  different  regions  and  time,  the  results  shown  LAGER  OPTICS  Vl.O  can 
automate  QC  process  under  various  conditions.  For  the  data  we  tested,  results  shown 
approximately  30.59%  of  the  total  profiles  were  sent  to  manual  QC  for  further  inspection,  when 
the  default  value  of  15%  flag  threshold  is  used.  A  tighter  threshold  will  prompt  more  data  to  be 
inspected  manually  (eg  10%  threshold  would  send  about  40%  data  for  manual  inspection,  based 
on  current  statistics  of  Slocum  data).  For  validation  tests,  a  default  value  of  10%  of  overall  flags, 
instead  of  15%  was  used,  to  ensure  a  tighter  than  usual  quality  control.  It  performs  extremely 
well  with  optics  data  from  slocum  glider  sll  13  during  Arabian  Gauntlet  2009  Exercise,  with 
under  1%  (actual  value  0.38%)  of  sampled  profiles  failed  to  be  checked  by  the  operator.  It  is 
worth  noticing  that  this  validation  test  was  done  using  more  questionable  profile  segments.  It 
should  be  anticipated  that  the  average  performance  should  be  better  than  shown  here,  both  in 
terms  of  erroneous  profiles  going  into  RTDHS  as  well  as  numbers  of  profiles  sent  to  MUG.  It 
has  been  observed  that  optical  data  from  Slocum  gliders  are  inherently  more  noisy  and  error- 
prone  for  various  reasons,  including  data  storage  and  transmission,  when  compared  to  those  from 
Seagliders.  Under  similar  conditions,  for  example,  less  than  1%  of  Seaglider  data  would  be  sent 
for  manual  inspection,  even  with  a  high,  10%  threshold.  (The  exact  figure  using  Table  5,  would 
be  0.1 1%  for  a  15%  threshold,  and  0.81%  for  10%  threshold.) 
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Table  1.  NRL/NAVO  glider  optical  parameters  for  BUFR  table. 

Traditional  optical  variable  names  and  units  are  used.  Typieal  values  listed  are  adapted  from  open  oeean  or  elear  water  environments. 
Bandwidth  parameters  are  from  current  or  past  instruments  tested  on  Sloeum  gliders  and  Seagliders,  with  names  denoted  from  raw 
data  stream.  The  speetral  information  is  not  listed  in  this  table  and  is  subjeet  to  ehange  depending  on  individual  sensors  onboard 
gliders.  Multiple  wavelengths  of  eertain  sensors  (e.g.  baekseattering  or  bb)  should  be  expeeted  and  have  been  implemented  in  the 
BUFR  table  by  means  of  repetition  (See  Table  2). 


Name 

Description 

Units 

Typical 

value 

Range 

Resolution 

(step) 

possible 
channel  #s 

Bandwidth 

(nm) 

Instruments 

bb 

baekseatter 

m-^ 

0.004 

0-0.1 

0.0001 

<255 

30 

bb2f,bb3slo 

Ed 

Irradianee 

wm'^ 

10 

0-500 

0.1 

<512 

10 

OCR504/7I 

(*) 

c 

beam 

attenuation 

m-^ 

0.2 

0-100 

0.01 

<255 

10 

BAM  (x) 

b 

Scattering 

m‘^ 

0.1 

0-100 

0.01 

<255 

10 

AUVb 

vis 

visibility 

m 

10 

0.1-100 

0.1 

<255 

- 

AUVb,SAM, 

BAM 

PAR 

photo,  avail,  rad. 

mEinstein. 

m-^ 

100 

0-2500 

1 

<255 

- 

ECO-PAR  (x) 

Fl-ehl 

[Chlorophyll-a] 

mg  m'"* 

0.5 

0-100 

0.01 

470/695 

10 

bb2f,fl3 

Fl-phyeo 

phycoerythrin 

ppb 

0.05 

0-100 

0.005 

540/570 

10 

fl3 

Fl-edom 

eolor  dissolved 

ppb  QS 

2 

0-2000 

0.1 

370/460 

10/70 

bb2f,fl3 
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Table  2.  Glider  optical  parameters  in  BUFR  format  for  LAGER  Optics  V1.0. 


F 

X 

Y 

Scale 

Reference 

Value 

Bit 

Width 

Unit 

Description 

0 

60 

001 

4 

0 

10 

m-1 

Backscatter 

0 

60 

002 

1 

0 

9 

wm-2 

Irradiance 

0 

60 

003 

2 

0 

14 

m-1 

Beam  Attenuation 

0 

60 

004 

2 

0 

14 

m-1 

Seattering 

0 

60 

005 

1 

0 

10 

m 

Visibility 

0 

60 

006 

0 

0 

12 

mEinstien.m-2 

Photosynthetic  Available  Radiation 

0 

60 

007 

2 

0 

14 

mg/m3 

Chlorophyll-a 

0 

60 

008 

3 

0 

17 

Ppb 

Cyanobaeteria 

0 

60 

009 

1 

0 

15 

Ppb 

Color  Dissolved 

0 

60 

010 

0 

0 

10 

Nm 

Center  Wavelength 

0 

60 

on 

0 

0 

7 

Nm 

Bandwidth 

0 

60 

012 

0 

0 

96 

CC1TT_1A5 

Instrument 
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Table  3.  BUFR  Quality  Control  (Class  63) 


F 

X 

Y 

Scale 

Reference 

Value 

Bit 

Width 

Unit 

Description 

0 

63 

001 

0 

0 

6 

Code  Table 

Automated  QC  Flags  -  Temperature 

0 

63 

002 

0 

0 

6 

Code  Table 

Automated  QC  Flags  -  Salinity 

0 

63 

003 

0 

0 

4 

Code  Table 

Manual  QC  Flags  -  Temperature 

0 

63 

004 

0 

0 

4 

Code  Table 

Manual  QC  Flags  -  Salinity 

0 

63 

005 

0 

0 

4 

Code  Table 

Full  Profile  QC  Flag  -  Temperature 

0 

63 

006 

0 

0 

4 

Code  Table 

Full  Profile  QC  Flag  -  Salinity 

0 

63 

007 

0 

0 

6 

Code  Table 

Automated  QC  Flags  -  Optics 

0 

63 

008 

0 

0 

4 

Code  Table 

Manual  QC  Flags  -  Optics 

0 

63 

009 

0 

0 

4 

Code  Table 

Full  Profile  QC  Flag  -  Backscatter 

0 

63 

010 

0 

0 

4 

Code  Table 

Full  Profile  QC  Flag  -  Irradiance 

0 

63 

on 

0 

0 

4 

Code  Table 

Full  Profile  QC  Flag  -  Beam  Attenuation 

0 

63 

012 

0 

0 

4 

Code  Table 

Full  Profile  QC  Flag  -  Scattering 

0 

63 

013 

0 

0 

4 

Code  Table 

Full  Profile  QC  Flag  -  Visibility 

0 

63 

014 

0 

0 

4 

Code  Table 

Full  Profile  QC  Flag  -  Photosynthetic  Available  Radiation 

0 

63 

015 

0 

0 

4 

Code  Table 

Full  Profile  QC  Flag  -  Chlorophyll-a 

0 

63 

016 

0 

0 

4 

Code  Table 

Full  Profile  QC  Flag  -  Cyanobacteria 

0 

63 

017 

0 

0 

4 

Code  Table 

Full  Profile  QC  Flag  -  Color  Dissolved 

0 

63 

018 

0 

0 

4 

Code  Table 

Full  Profile  QC  Flag  -  Center  Wavelength 
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Table  4.  Statistical  Results  under  Different  Testing  Conditions 


Data  from  Slocum  gliders  are  used  to  test  the  performanee  of  current  Glider  Optics  Vl.O  auto  QC  algorithms.  Notice  that  optieal  data 
from  Slocum  gliders  typically  suffer  from  more  noises  and  data  issues  that  those  from  Seagliders.  This  is  the  main  reason  for  choosing 
Slocum.  Detailed  eomparisons  between  optics  from  Slocum  and  Seagliders  can  be  found  from  Table  5  and  more  details  in  text. 


Glider 

Date 

Exercise 

Number  of 
Profiles 

30% 

flagged 

25% 

flagged 

20% 

flagged 

15% 

flagged 

10% 

flagged 

5% 

flagged 

S1113 

Aug  2009 

Arabian  Gauntlet 

1038 

0.289 

0.624 

2.70 

12.1 

40.4 

83.8 

S1005 

June  2008 

BioSpace 

360 

3.06 

3.61 

6.67 

19.4 

37.22 

80.3 

S1016 

June  2008 

BioSpaee 

168 

0 

0 

0 

0.595 

0.595 

0.595 

S1027 

June  2008 

BioSpace 

216 

63.0 

63.0 

63.0 

63.0 

67.1 

67.1 

S1079 

June  2008 

BioSpace 

1188 

37.1 

37.1 

37.1 

37.1 

40.6 

71.5 

S1082 

June  2008 

BioSpace 

1454 

0.138 

0.138 

1.03 

2.82 

11.9 

61.9 

S183 

June  2008 

BioSpace 

408 

60.3 

60.3 

60.3 

60.3 

60.3 

80.6 

S1016 

July  2008 

RIMPAC 

2235 

47.7 

47.7 

48.2 

49.4 

50.9 

54.5 

Total 

7067 

26.45 

26.56 

27.38 

30.59 

38.63 

62.54 
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Table  5.  Overall  flag  stats  for  all  profiles  tested 


lArea  ' 

Glider 

Threshold 

NoFiles(%) 

Manual 

QC(%) 

Individual  Flag  % 

2 

4 

5 

6 

7 

GL_PGulf_200908 

SI113 

1038 

30 

0.29 

GL  PGulf  200908 

SI113 

1038 

25 

0.67 

GL_PGulf_200908 

SI113 

1038 

20 

2.7 

GL  PGulf  200908 

SI113 

1038 

15 

12.14 

0.273 

5.28 

1.451 

0.013 

0.26 

GL_PGulf_200908 

SI113 

1038 

10 

40.37 

GL_PGulf_200908 

SI113 

1038 

5 

83.82 

GL_BIOSPACE_200806 

SI005 

168 

30 

0 

GL_BIOSPACE_200806 

SI005 

168 

25 

0 

GL  BIOSPACE  200806 

SI005 

168 

20 

0 

GL_BIOSPACE_200806 

SI005 

168 

15 

0.6 

0.034 

3.097 

1.929 

2.323 

1.069 

GL_BIOSPACE_200806 

slOOS 

168 

10 

0.6 

GL_BIOSPACE_200806 

slOOS 

168 

5 

19.64 

GL_BIOSPACE_200806 

SI005 

168 

2 

82.74 

GL_BIOSPACE_200806 

SI016 

360 

30 

3.06 

GL_BIOSPACE_200806 

SI016 

360 

25 

3.61 

GL  BIOSPACE  200806 

SI016 

360 

20 

6.67 

GL_BIOSPACE_200806 

SI016 

360 

15 

19.44 

0.107 

0.601 

1.621 

0 

0.07 

GL  BIOSPACE  200806 

SI016 

360 

10 

37.22 

GL_BIOSPACE_200806 

SI016 

360 

5 

80.28 

GL_BIOSPACE_200806 

SI077 

219 

30 

63.01 

GL  BIOSPACE  200806 

SI077 

219 

25 

63.01 

GL_BIOSPACE_200806 

SI077 

219 

20 

63.01 

GL  BIOSPACE  200806 

SI077 

219 

15 

63.01 

0.037 

1.061 

0.978 

24.404 

0.627 

GL_BIOSPACE_200806 

SI077 

219 

10 

67.12 

GL_BIOSPACE_200806 

SI077 

219 

5 

79 

GL  BIOSPACE  200806 

SI079 

1188 

30 

37.12 

GL_BIOSPACE_200806 

SI079 

1188 

25 

37.12 

GL  BIOSPACE  200806 

SI079 

1188 

20 

37.12 

GL_BIOSPACE_200806 

SI079 

1188 

15 

37.12 

0.005 

2.375 

1.065 

13.546 

0.615 

GL  BIOSPACE  200806 

SI079 

1188 

10 

40.57 

GL_BIOSPACE_200806 

SI079 

1188 

5 

71.55 

GL_BIOSPACE_200806 

SI082 

1454 

30 

0.14 

GL  BIOSPACE  200806 

SI082 

1454 

25 

0.14 

GL  BIOSPACE  200806 

SI082 

1454 

20 

1.03 

24 


GL_ 

CD 

O 

00 

O 

O 

LU 

o 

< 

CL 

CO 

O 

ca 

SI082 

1454 

15 

GL_ 

_BIOSPACE_200806 

SI082 

1454 

10 

GL_ 

_BIOSPACE_200806 

SI082 

1454 

5 

GL_ 

CD 

O 

00 

O 

O 

LU 

O 

< 

CL 

CO 

O 

CO 

SI083 

408 

30 

GL_ 

_BIOSPACE_200806 

SI083 

408 

25 

GL_ 

_BIOSPACE_200806 

SI083 

408 

20 

GL_ 

_BIOSPACE_200806 

SI083 

408 

15 

GL_ 

_BIOSPACE_200806 

SI083 

408 

10 
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Appendix  A 

opticsvariablesinfo .  dat 


This  file  is  the  parameter  file  for  most  of  the  automated  opties  QC  modules,  as  shown  in  Figure 
2.  This  file  resides  in  SETUP  FILES  subdireetory  under  EAGER.  Modifieations  of  parameters 
in  this  file  will  directly  affect  the  program  performance.  It  should  be  handled  with  caution,  and 
knowledge  of  specific  optical  sensors  and  its  calibration  procedures. 


# 

#  File  Format  Description 

# 

#  General 

#  All  lines  must  start  in  column  1. 

# 

#  Lines  starting  with  asterisk  (*)  are  comments 

#  Comment  lines  can  occur  anywhere  within  the  file  and  are  not 

#  required  to  occur  anywhere  within  the  file.  For  readability,  a 

#  comment  line  should  be  placed  between  each  block  of  lines  defining 

#  a  variable  or  variable  group. 

# 

#  Lines  with  a  variable  definition  block  can  be  in  any  order. 

# 

#  If  a  block  of  lines  following  the  line:  groupsize  N,  where  N 

#  is  greater  than  1,  must  have  N  variable,  type,  name,  units, 

#  and  validrange  lines.  In  addition,  if  there  is  at  least  one 

lines . 


# 

#  Description  of 

# 

#  Line  Type 

#  1 .  nvarinput 

# 

# 

# 

# 

line . 

# 

#  2 .  nvaroutput 

# 

# 

#  3 .  inname 

# 

# 

# 

# 

#  4.  calibrated 

# 

# 

# 

# 

#  5 .  type 


each  line  type 
Description 

Number  of  variables  input  as  processing  group. 
Example:  groupinsize  1 

groupsize  followed  by  integer  indicating 

the  number  of  variable  required  during  processing. 

Each  new  group  sequence  MUST  begin  with  a  groupinsize 


Number  of  variables  output  from  processing  group. 
Example:  groupoutsize  2 

If  missing,  then  groupoutsize  I  assumed, 
input  variable  name 
Example : 

If  nvarinput  >  1,  then  each  variable  name  in  group  must 
be  listed  on  separate  lines. 

Variable  is  already  calibrated. 

Examples : 
calibrated  yes 
calibrated  no 

Generic  variable  type,  e.g.,  bb,  ed,  FL-chl 
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Units  of  calibrated  variable. 


# 

#  6.  units 

# 

#  7.  outname  Name  of  calibrated  variable. 

# 

#  8.  validrangemin  Valid  range  minimum  of  calibrated  variable 

# 

#  9.  validrangemax  Valid  range  maximum  of  calibrated  variable 

# 

#  9.  wavelength 

# 

#  10.  bandwidth 

# 

#  11.  scalefactor 

# 

#  12.  enddef  Each  group  sequence  MUST  end  with  an  enddef  line. 

# 

nvarinput  =  1 
inname  =  wlbb2f_blueCount 
calibrated  =  no 
instrument  =  bb2f 

description  =  backscattering  coefficient 
type  =  bb 
units  =  1/m 
outname  =  backscatter 
validrangemin  =  0.0 
validrangemax  =  0.1 
resolution  =  0.0001 
wavelength  =  470 
bandwidth  =  30 
scalefactor  =  1 
spike  =  l.e-3 
stdevfactor  =  2.0 
abslimfactor  =  3.e-4 
enddef 
# 

nvarinput  =  1 
inname  =  wlbb2f  redCount 
calibrated  =  no 
instrument  =  bb2f 

description  =  backscattering  coefficient 
type  =  bb 
units  =  1/m 
outname  =  backscatter 
validrangemin  =  0.0 
validrangemax  =  0.1 
resolution  =  0.0001 
wavelength  =  700 
bandwidth  =  30 
scalefactor  =  1 
spike  =  l.e-3 
stdevfactor  =  2.0 
abslimfactor  =  3.e-4 
enddef 
# 

nvarinput  =  1 

inname  =  wlbb2f  fluorCount 
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calibrated  =  no 
instrument  =  bb2f 

description  =  chlorophyll-a  concentration  from  fluorescence 
type  =  Flchl 
units  =  mg/m**3 
outname  =  chlorophyll-a 
validrangemin  =  0 
validrangemax  =  100 
resolution  =  0.01 
wavelength  =  470 
scalefactor  =  1 
spike  =  0.2 
stdevfactor  =  2.0 
abslimfactor  =  0.03 
enddef 
# 

nvarinput  =  1 

inname  =  sci  bb3slo  b470  scaled 
calibrated  =  yes 
instrument  =  bb3slo 

description  =  backscattering  coefficient 
type  =  bb 
units  =  1/m 
outname  =  backscatter 
validrangemin  =  0 
validrangemax  =  0.1 
resolution  =  0.0001 
wavelength  =  470 
scalefactor  =  1 
spike  =  l.e-3 
stdevfactor  =  2.0 
abslimfactor  =  3.e-4 
enddef 
# 

nvarinput  =  1 

inname  =  sci  bb3slo  b532  scaled 
calibrated  =  yes 
instrument  =  bb3slo 

description  =  backscattering  coefficient 
type  =  bb 
units  =  1/m 
outname  =  backscatter 
validrangemin  =  0 
validrangemax  =  0.1 
resolution  =  0.0001 
wavelength  =  532 
scalefactor  =  1 
spike  =  l.e-3 
stdevfactor  =  2.0 
abslimfactor  =  3.e-4 
enddef 
# 

nvarinput  =  1 

inname  =  sci_bb3slo_b660  scaled 
calibrated  =  yes 
instrument  =  bb3slo 

description  =  backscattering  coefficient 
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type  =  bb 
units  =  l/m 
outname  =  backscatter 
validrangemin  =  0 
validrangemax  =  0.1 
resolution  =  0.0001 
wavelength  =  660 
scalefactor  =  1 
spike  =  l.e-3 
stdevfactor  =  2.0 
abslimf actor  =  3.e-4 
enddef 
# 

nvarinput  =  1 

inname  =  sci  ocr504I  irradl 
calibrated  =  yes 
instrument  =  ocr504I 
description  =  downward  irradiance 
type  =  ed 
units  =  w/m**2 
outname  =  irradiance 
validrangemin  =  0 
validrangemax  =  500 
resolution  =  0.1 
wavelength  =  0 
bandwidth  =  10 
scalefactor  =  1.0 
stdevfactor  =  3.0 
abslimfactor  =  0.3 
spike  =  50 
enddef 
# 

nvarinput  =  1 

inname  =  sci_ocr504I  irrad2 
calibrated  =  yes 
instrument  =  ocr504I 
description  =  downward  irradiance 
type  =  ed 
units  =  w/m**2 
outname  =  irradiance 
validrangemin  =  0 
validrangemax  =  500 
resolution  =  0.1 
wavelength  =  0 
bandwidth  =  10 
scalefactor  =  1.0 
stdevfactor  =  3.0 
abslimfactor  =  0.3 
spike  =  50 
enddef 
# 

nvarinput  =  1 

inname  =  sci  ocr504I  irrad3 
calibrated  =  yes 
instrument  =  ocr504I 
description  =  downward  irradiance 
type  =  ed 
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units  =  w/m**2 
outname  =  irradiance 
validrangemin  =  0 
validrangemax  =  500 
resolution  =  0.1 
wavelength  =  0 
bandwidth  =  10 
scalefactor  =  1.0 
stdevfactor  =  3.0 
abslimf actor  =  0.3 
spike  =  50 
enddef 
# 

nvarinput  =  1 

inname  =  sci  ocr504I  irrad4 
calibrated  =  yes 
instrument  =  ocr504I 
description  =  downward  irradiance 
type  =  ed 
units  =  w/m**2 
outname  =  irradiance 
validrangemin  =  0 
validrangemax  =  500 
resolution  =  0.1 
wavelength  =  0 
bandwidth  =  10 
scalefactor  =  1.0 
stdevfactor  =  3.0 
abslimfactor  =  0.3 
spike  =  50 
enddef 
# 

nvarinput  =  1 

inname  =  sci  bbfl2s_bb  sig 
calibrated  =  no 
instrument  =  bbfl2s 

description  =  backscattering  coefficient 
type  =  bb 
units  =  1/m 
outname  =  backscatter 
validrangemin  =  0 
validrangemax  =  0.0001 
resolution  =  0.0001 
wavelength  =  0 
bandwidth  =  30 
scalefactor  =  1 
spike  =  l.e-3 
stdevfactor  =  2.0 
abslimfactor  =  0.3 
enddef 
# 

nvarinput  =  1 

inname  =  sci  bbfl2s  chlor  sig 
calibrated  =  no 
instrument  =  bbfl2s 

description  =  chlorophyll-a  concentration  from  fluorescence 
type  =  Flchl 
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units  =  mg/m**3 
outname  =  chlorophyll-a 
validrangemin  =  0 
validrangemax  =  100 
resolution  =  0.01 
wavelength  =  0 
bandwidth  =  10 
scalefactor  =  1 
spike  =  0.1 
stdevfactor  =  2.0 
abslimf actor  =  0.03 
enddef 
# 

nvarinput  =  1 

inname  =  sci  bbfl2s  cdom  sig 
calibrated  =  no 
instrument  =  bbfl2s 

description  =  color  dissolved  organic  matter  fluorescence 
type  =  Flcdom 
units  =  ppb  QS 
outname  =  color  dissolved 
validrangemin  =  0 
validrangemax  =  2000 
resolution  =  0.1 
wavelength  =  0 
bandwidth  =  10 
scalefactor  =  1 
spike  =  5.0 
stdevfactor  =  2.0 
abslimfactor  =  0.3 
enddef 
# 

nvarinput  =  1 

inname  =  sci  bbfl2s_bb  scaled 
calibrated  =  yes 
instrument  =  bbfl2s 

description  =  backscattering  coefficient 
type  =  bb 
units  =  1/m 
outname  =  backscatter 
validrangemin  =  0 
validrangemax  =  0.1 
resolution  =  0.0001 
wavelength  =  661 
bandwidth  =  30 
scalefactor  =  1 
spike  =  l.e-3 
stdevfactor  =  2.0 
abslimfactor  =  3.e-4 
enddef 
# 

nvarinput  =  1 

inname  =  sci  bbfl2s  chlor  scaled 
calibrated  =  yes 
instrument  =  bbfl2s 

description  =  chlorophyll-a  concentration  from  fluorescence 
type  =  Flchl 
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units  =  mg/m**3 
outname  =  chlorophyll-a 
validrangemin  =  0 
validrangemax  =  100 
bandwidth  =  30 
wavelength  =  0 
resolution  =  0.01 
scalefactor  =  1 
spike  =  0.1 
stdevfactor  =  2.0 
abslimf actor  =  0.03 
enddef 
# 

nvarinput  =  1 

inname  =  sci  bbfl2s  cdom  scaled 
calibrated  =  yes 
instrument  =  bbfl2s 

description  =  color  dissolved  organic  matter  fluorescence 
type  =  Flcdom 
units  =  ppb  QS 
outname  =  color  dissolved 
validrangemin  =  0 
validrangemax  =  2000 
bandwidth  =  10 
resolution  =  0.1 
wavelength  =  0 
scalefactor  =  1 
spike  =  1.0 
stdevfactor  =  2.0 
abslimfactor  =  0.3 
enddef 
# 

nvarinput  =  2 
nvaroutput  =  3 
inname  =  sci  auvb  ref 
inname  =  sci  auvb  sig 
calibrated  =  no 
instrument  =  auvb 
instrument  =  auvb 
instrument  =  auvb 

description  =  derived  diver  visibility  range 

description  =  beam  attenuation  coefficient 

description  =  total  scattering  coefficient 

type  =  vis 

type  =  c 

type  =  b 

units  =  m 

units  =  1/m 

units  =  1/m 

outname  =  visibility 

outname  =  attenuation 

outname  =  scattering 

scalefactor  =  1 

scalefactor  =  1 

scalefactor  =  1 

wavelength  =  0 

wavelength  =  0 
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wavelength  =  0 
resolution  =  0.01 
resolution  =  0.01 
resolution  =  0.01 
validrangemin  =  0.1 
validrangemax  =  100 
validrangemin  =  0 
validrangemax  =  100 
validrangemin  =  0 
validrangemax  =  100 
bandwidth  =  10 
bandwidth  =  10 
bandwidth  =  10 
spike  =  1.0 
spike  =  1.0 
spike  =  1.0 
stdevfactor  =  2.0 
stdevfactor  =  2.0 
stdevfactor  =  2.0 
abslimfactor  =  0.03 
abslimfactor  =  0.03 
abslimfactor  =  0.03 
enddef 
# 

nvarinput  =  1 

inname  =  sci  fl3slo  chlor  units 
calibrated  =  yes 
instrument  =  fl3slo 

description  =  chlorophyll-a  concentration  from  fluorescence 
type  =  Flchl 
units  =  mg/m**3 
outname  =  chlorophyll-a 
scalefactor  =  1 
resolution  =  0.01 
wavelength  =  0 
bandwidth  =  30 
validrangemin  =  0 
validrangemax  =  100 
spike  =  0.1 
stdevfactor  =  2.0 
abslimfactor  =  0.03 
enddef 
# 

nvarinput  =  1 

inname  =  sci  fl3slo  phyco  units 
calibrated  =  yes 
instrument  =  fl3slo 

description  =  phycoerytherin  fluorescence 

type  =  Flphyco 

units  =  ppb 

scalefactor  =  1 

outname  =  cyanobacteria 

resolution  =  0.03 

wavelength  =  0 

bandwidth  =  10 

validrangemin  =  0 

validrangemax  =  2000 


36 


spike  =  0.3 
stdevfactor  =  2.0 
abslimf actor  =  0.09 
enddef 
# 

nvarinput  =  1 

inname  =  sci  fl3slo  cdom  units 
calibrated  =  yes 
instrument  =  fl3slo 

description  =  color  dissolved  organic  matter  fluorescence 

type  =  Flcdom 

units  =  ppb  QS 

outname  =  color  dissolved 

scalefactor  =  1 

resolution  =  0.09 

wavelength  =  0 

bandwidth  =  10 

validrangemin  =  0 

validrangemax  =  100 

spike  =  0.9 

stdevfactor  =  2.0 

abslimfactor  =  0.27 

enddef 
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